On the definition of temperature in dense granular media 
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In this Letter we report the measurement of a pseudo-temperature for compacting granular media 
on the basis of the Fluctuation-Dissipation relations in the aging dynamics of a model system. From 
the violation of the Fluctuation-Dissipation Theorem an effective temperature emerges (a dynamical 
temperature Tdyn) whose ratio with the equilibrium temperature TJ' depends on the particle density. 
We compare the results for the Fluctuation-Dissipation Ratio (FDR) Tdyn/T^'' at several densities 
with the outcomes of Edwards' approach at the corresponding densities. It turns out that the FDR 
and the so-called Edwards' ratio coincide at several densities (very different ages of the system), 
opening in this way the door to experimental checks as well as theoretical constructions. (PACS: 
05.70.Ln, 05.20.-y, 45.70.Cc) 
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The study of compact granular matter through statisti- 
cal physics tools is the subject of a sustained interest [Q . 
Granular media enter only partially in the framework 
of equilibrium statistical mechanics and their dynamics 
constitutes a very complex problem of non-equilibrium, 
which poses novel questions and challenges to theorists 
and experimentalists. The very possibility to construct 
a coherent statistical mechanics for these systems is still 
matter of debate, although everybody agrees that such 
an approach, if possible, would allow for a much deeper 
and global understanding of the problem. 

One of the main obstacles in this direction is the non- 
thermal character of these systems: thermal energy is 
so negligibly small with respect to other energy contri- 
butions (e.g. potential energy) that for all the practical 
purposes these systems live virtually at zero tempera- 
ture. One of the most important consequences is that, 
unless perturbed in some way (e.g. driving energy into 
the system), a granular system cannot explore sponta- 
neously its phase space but it remains trapped in one of 
the numerous metastable configurations. Understanding 
the structure of the phase space which is left invariant 
by the dynamics is then crucial for the construction of 
a thermodynamical description of these non-thermal sys- 
tems. 

A very ambitious approach, in this direction, has been 
put forward by S. Edwards and co-workers |§, D, by 
proposing an equivalent of the microcanonical ensemble: 
macroscopic quantities in a jammed situation should be 
obtained by a flat average over all blocked configurations 
(i.e. in which every grain is unable to move) of given vol- 
ume, energy, etc... The strong assumption here is that 
all blocked configurations are treated as equivalent and 
have the same weight in the measure. This approach, 
based on the idea of describing granular material with 
a small number of parameters, leads to the introduction 
of an entropy Sedw , given by the logarithm of the num- 
ber of blocked configurations of given volume, energy. 



etc., and its corresponding density Sedw = Sedw/N. As- 
sociated with this entropy are the state variables such 
as 'compactivity' — -^SedwiV) and 'temperature' 

T ^ — -^SedwiE). 

Very recently, important progresses in this direction 
have been reported in various contexts: a tool to system- 
atically construct Edwards' measure, defined as the set 
of blocked configurations of a given model, was proposed 
in 1^, ^ ; it was used to show that the outcome of the ag- 
ing dynamics of the Kob- Andersen model (a kinetically 
constrained lattice gas model) was correctly predicted by 
Edwards' measure. Moreover, the validity and relevance 
of Edwards' measure have been demonstrated for one- 
dimensional phenomenological models for spin mod- 
els with "tapping" dynamics 0|, and for sheared hard 
spheres 

In this Letter we focus on the definition of a 
pseudo-temperature for granular media on the basis 
of the Fluctuation-Dissipation relations in the out-of- 
cquilibrium, aging, dynamics of a model system and 
on its relation to Edwards' measure. From the violation 
of the Fluctuation-Dissipation Theorem an effective tem- 
perature emerges (from now onward indicated as dynam- 
ical temperature Tdyn) whose ratio with the equilibrium 
temperature T^'^ depends on the particle density. We 
compare the results for the Fluctuation-Dissipation Ratio 
(FDR) Tdyn I T^'^ at several densities with the outcomes 
of Edwards' approach at the corresponding densities. It 
turns out that the FDR and the so-called Edwards' ratio 
coincide at several densities (very different ages of the 
system), opening in this way the door to experimental 
checks as well as theoretical constructions. 

It is interesting to mention recent approaches that are 
complementary to ours. On the one hand, in the context 
of aging supercooled liquids, the inherent structure strat- 
egy does not address the question of a determination of a 
static distribution, but promising results show that this 
measure, if it exists, is insensitive to the details of the 
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thermal history . The link between this strategy and 
Edwards' measure has been discussed in On the 

other hand, recent works on the possibility of a dynam- 
ical definition of temperature have focused on sheared, 
stationary systems These studies are clearly 

complementary to the present one, which addresses the 
problem of the relation between the value of Tdyn and a 
static measure, in aging (non-stationary) systems. 

The model we consider is a version of the "Tetris" 
model H] , which has been shown to reproduce several fea- 
tures of granular media like aging |l^ , memory p6[ , 
self-structuring etc. In the framework of this model, 
some of us have already provided one of the first evi- 
dences of the validity of Edwards' measure ||] . 

We focus in particular on an homogeneous system with 
no preferential direction in order to avoid a ny kind of in- 
stability or large-scale structure formation . The case 
with gravity which imposes a preferential direction will 
be discussed elsewhere . In the version of the model 
we use, "T" -shaped particles diffuse on a square lattice, 
with the only constraint that no superposition is allowed: 
for two nearest-neighbor particles, the sum of the arms 
oriented along the bond connecting the two particles has 
to be smaller than the bond length (for each particle, 
the three arms of the 'T" have length jd, where d sets 
the bond size on the square lattice). The maximum den- 
sity allowed is then Pmax = 2/3. This model represents a 
clear example of a non-thermal system. The Hamiltonian 
is zero and the temperature itself is therefore irrelevant 
at equilibrium, only its ratio with an imposed chemical 
potential being important. 

The out-of-equilibrium compaction dynamics without 
gravity is implemented as follows: starting from an empty 
lattice, particles are randomly deposited, without diffu- 
sion and without violation of the geometrical constraints. 
This random sequential absorption process yields a re- 
producible initial density of p ~ 0.547. Alternating dif- 
fusions and additions of particles are then attempted, 
allowing to increase the density of the system, which re- 
mains homogeneous in the process. In order to overcome 
the problem related to the simulation of slow processes 
and obtain a reasonable number of different realizations 
to produce clean data, we have devised a fast algorithm 



(in the spirit of Bortz-Kalos-Lebowitz algorithm |19) 



where the essential ingredient is the updating of a list of 
mobile particles (whose number is rimob)- At each time 
step one selects a mobile particle and move it in the direc- 
tion chosen only if no violation of the constraints occurs. 
The time is incremented of an amount At = l/rimob- We 
refer to for the details of the fast algorithm. 

We have simulated lattices of linear size L = 
50, 100, 200, in order to ensure that finite-size effects 
were irrelevant. We have chosen periodic boundary con- 
ditions on the lattice, having verified that other types 
of boundary conditions (e.g closed ones) gave the same 
results. 



During the compaction, we monitor the following 
quantities: the density of particles p{t), the density of 
mobile particles Pmob{t)- Moreover in order to estab- 
lish the Fluctuation-Dissipation relations we measure the 
mean square displacement B{t + tw,tw) and the inte- 
grated response function xit + tw,tw)- The mean square 
displacement is defined as 

1 ^ 
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where N is the number particles present in the system 
at time tw, ri is the coordinate {x or y) of the i-th parti- 
cle and the brackets () indicate the average over several 
realizations. 

In order to measure the integrated response function 
xit tw,tw), we make a copy of the system at time 
and apply to it a small random perturbation, varying the 
diffusion probability of each particle from p — j to p"^ — 
j + f[ ■ e, where = ±1 is a random variable associated 
to each grain independently for each possible direction 
(r = x,y), and e represents the perturbation strength. 
For a constant field we obtain (see also Gl pi p8|): 
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where Anit + t^) = [t + ty,) - r^it + t^) is the 
difference between the displacements taking place in the 
two systems evolving with the same succession of random 
numbers. In this letter we present the results obtained 
with a perturbation strength e = 0.005, having checked 
that for 0.002 < e < 0.01 non-linear effects are absent. 

If the system was in equilibrium we would expect B 
and X to be linearly related by 



X{t ~^ tuntw) — rpeq B(t -\- tm,t-i^"), 



(1) 



where Xdyn is the so-called Fluctuation-Dissipation Ra- 
tio (FDR) which is unitary in equilibrium. Any de- 
viations from this linear law signals a violation of the 
Fluctuation-Dissipation Theorem (FDT) . Nevertheless it 
has been shown, first in mean- field models then 
in various simulations ||2l] , how in several aging sys- 
tems violations from (|l|) reduce to the occurrence of two 
regimes: a quasi-equilibrium regime with Xdyn = 1 (and 
time-translation invariance) for "short" time separations 
{t ^ tw), and the aging regime with Xdyn < 1 for large 
time separations. This second slope is typically referred 
to as a dynamical temperature Tdyn > Tj"^ such that 

Xdyn — T^'^ /Tdyn [ p3[ . 

In our case, as the density increases, an aging be- 
haviour is obtained: the system falls out of equilibrium 
and pmob{t) gets smaller than the corresponding value 
at equilibrium Accordingly, violations of (^ are ex- 
pected. 
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If the compaction process is stopped at a certain time 
tw, the system relaxes toward equilibrium and one ob- 
tains a time-translation invariant behaviour for x and B; 
this is the so-called regime of interrupted aging which 
features an increase of Pmob to its equilibrium value and 
a single linear relation for the x vs. B parametric plot. 
These measures therefore give us the value of the equi- 
librium fluctuation-dissipation ratio, T^'' , which does ac- 
tually not depend on the density reached at tw 

If, on the other hand, the compaction process is not 
stopped, the system features an aging behaviour and the 
X vs. B parametric plots, displayed in figure 0, show two 
different linear behaviours: after a first quasi-equilibrium 
regime in which the plot has slope T^"^, a violation of 
FDT is observed, with the existence of a dynamical tem- 
perature Tdyn which depends on and therefore on the 
density. We obtain in particular X^j,„ = 0.646 ± 0.002, 
^lyn = 0.767 ±0.005, Xjy^ = 0.784 ±0.005 at tl = 10^, 
^2 ^ 3 • 10^ and tl=5- 10^. 



become blocked and therefore equilibrium and Edwards' 
measures become equivalent. 
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FIG. 2: Edwards' ratio X^dw as a function of density. The 
horizontal lines correspond to the dynamical ratios Xdyn, 
measured (from bottom to top) at t^j = 10*, 3 ■ 10*, 5 • 10* 
and determine the values pi ~ 0.596, p2 ~ 0.603, pa ~ 0.605, 
to be used in the comparison with the results reported in Fig. 
3. 
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FIG. 1: Einstein relation in the Tetris model: plot of the 
mobility xi^w + 't,tw)T^'' vs. the mean-square displacement 
B{tw+t, tw), for — 10* and = 3 - 10*. The unitary slope 
of the full straight line corresponds to the equilibrium case, 
obtained for the dynamics at constant density (interrupted 
aging). 



We are now able to compare the values of the dynam- 
ical measures with the outcome of Edwards' measure. 
The construction of the equilibrium and Edwards' mea- 
sures has been described in Q. In particular, an efficient 
sampling of the blocked configurations, and therefore Ed- 
wards' measure, is obtained by the use of an auxiliary 
model whose energy is defined to be the number of mo- 
bile particles: the introduction of an auxiliary tempera- 
ture and an annealing procedure then yields configura- 
tions with no mobile particles. Having computed equi- 
librium and Edwards' entropies as a function of density, 
as in @, we measure the ratio of the slopes (that we 
denote as Edwards' ratio) 



XEdw — 



dsEdwjp) / dSequil (fi) 



dp 



dp 



(2) 



which is plotted in figure (g). This ratio approaches 1 as 
p — > 2/3, since at the maximum density all configurations 



The values of the dynamical Fluctuation-Dissipation 
ratio are also reported in the same figure (with no error 
bars since they are too small) and yield the following 
densities: pi « 0.596 for t]^ = 10*, p2 « 0.603 for = 
3.10", p3 « 0.605 for tl, = 10^. On the other hand, 
the evolution of the density of the system during the 
measurements is reported in figure (H). 







density 


density interval 


= 10* 


0.646 


0.596 


[0.584, 0.597] 


t» = 3 ■ 10* 


0.767 


0.603 


[0.599, 0.605] 


t„ = 5 • 10* 


0.784 


0.605 


[0.603, 0.606] 



TABLE I; Fluctuation-Dissipation Ratio {Xdyn) as obtained 
from numerical data fits with the corresponding densities ob- 
tained from figure (^). The last column reports the density 
intervals explored during the compaction dynamics (starting 
from the times where deviations from equilibrium become ev- 
ident) for U = 10*, 3 • 10*, 5 • 10*. 



Since the measurements are performed during the com- 
paction, the density is evolving, going from p{t^) to 
p{tw + tmax)- In each case, we obtain that indeed 
Pi G [p(^L))7 + tmax)] where we have denoted with 
Pi the densities obtained from figure ^ for different val- 
ues of {i = 1,2,3). More precisely, pi is very close 
to p(t^ + tmax)- This is to be expected since the mea- 
sure of the FDT violation is made for times much larger 
than and, since the compaction is logarithmic, the 
system spends actually more time at densities close to 
p(C + tmax) than to p(i^). 

The validation of Edwards' hypothesis in various model 
systems has made important steps forwards recently; 
here we have focused, for a model with only geometri- 
cal constraints, on the definition of a dynamical temper- 
ature and on its link with Edwards' measure. While the 
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FIG. 3: Evolution of the density during the measurements 
of X and B, for t„ = 10"*, 3 • lO*, 5 ■ lO*. The evolution 
during the quasi-equilibrium part is plotted with lines, and 
during the violation of FDT with symbols. The horizontal 
lines correspond to the densities pi, p2, ps from Fig 2. 



density increases, the measured dynamical temperature 
decreases, following closely the ratio between the equi- 
librium and the Edwards' entropies, the latter being ob- 
tained through a flat sampling of blocked configurations. 
While Edwards' proposal is supposed to be only valid 
asymptotically, i.e. when one-time quantities are almost 
stationary, our study clearly shows that it actually yields 
good results even in a pre-asymptotic regime, when the 
density is still evolving a lot. 



T NO remarks arc in order. It is interesting to invcsti- 
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